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THE FOKKER-PLANCK EQUATION FOR BOSONS IN 2D: 
WELL-POSEDNESS AND ASYMPTOTIC BEHAVIOR 


JOSE A. CANIZO, JOSE A. CARRILLO, PHILIPPE LAURENgOT, AND JESUS ROSADO 


Abstract. We show that solutions of the 2D Fokker-Planck equation for bosons 
are defined globally in time and converge to equilibrium, and this convergence is 
shown to be exponential for radially symmetric solutions. The main observation 
is that a variant of the Hopf-Cole transformation relates the 2D equation in 
radial coordinates to the usual linear Fokker-Planck equation. Hence, radially 
symmetric solutions can be computed analytically, and our results for general 
(non radially symmetric) solutions follow from comparison and entropy argu¬ 
ments. In order to show convergence to equilibrium we also prove a version 
of the Csiszar-Kullback inequality for the Bose-Einstein-Fokker-Planck entropy 
functional. 
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1. Introduction 


Introduced in the early 1900s [28, 16], Fokker-Planck equations have demon¬ 
strated to be a powerful tool to study macroscopic qualitative properties of systems 
of many interacting particles. They are particularly relevant in the held of sta¬ 
tistical mechanics since the probability laws of Ornstein-Uhlenbeck or Langevin 
processes follow this type of partial differential equation [21, 29, 19, 17]. They 
are usually obtained as grazing collision limits of the Boltzmann equation or ex¬ 
pansions from the master equation in classical [23, 12] and in quantum contexts 
[22, 24, 27, 25, 30, 2], and hnd applications in multiple areas [10, 6, 20]. In 
this way, nonlinear Fokker-Planck equations appear naturally as simplihed ki¬ 
netic equations with typically the same equilibrium states as the more physically 
grounded Boltzmann-type kinetic or master equations, but keeping some of the 
essential features and difficulties. 
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In this paper, we work on the particular case of the model introduced in the 
nineties by Kaniadakis and Quarati [26] to describe the distribution of quantum 
particles in a spatially homogeneous setting leading to possible concentrations in 
hnite time. From the physics point of view, it is quite relevant due to its relation to 
the Bose-Einstein condensate phenomena [14, 15, 31]. Different models have also 
been proposed to study this process, mostly based upon the Schrodinger equation 
and its Gross-Pitaevskii variants [1, 13, 32]. 

The basic idea is to introduce a nonlinear Fokker-Planck operator such that the 
Bose-Einstein distributions are their only possible steady states. These distribu¬ 
tions are the family of functions of the form 



1 

/3exp(|n|72) - 1’ 


V e 


( 1 ) 


where /3 > 1 is a parameter. They obviously satisfy the equation 


V„log(D7')(„) + „ = 0 ^ VJ>)+vfi(l + f!^)(v) = 0 

V 1 T Joo / 

for all /5 > 1. Clearly, for each (3 > 1 is a stationary solution to the following 
initial value problem for / = f{t,v), with f > 0 and v G 

dtf = Af+ div{vf{l +f)), t>0, neM'^, (2a) 

f{0,v) = fo{v), ugM'^. (2b) 

The mass at time f of a solution / to eq. (2) is dehned by its integral, and it is a 
conserved quantity of the time evolution, that is. 


f(t, v) dv = m 


fo{v) dn , f > 0 . 


Problem (2) is referred to as the Cauchy problem for the Bose-Einstein-Fokker- 
Planck (BEFP) equation or the Fokker-Planck equation for bosons. It is easy to 
check that the functions with (3 > 1 given by (1) are the only possible finite-mass 
smooth stationary solutions of (2). 

In dimension d > 3 there are no stationary solutions with mass larger than that 
of while in dimensions d = 1,2 the mass of is infinite and there is exactly 
one stationary solution for each fixed hnite mass. We dehne the critical mass me 
by 

which is hnite in dimensions d> 3 and equal to cxd if d = 1, 2. The mass of is 
actually explicitly computable in dimension d = 2: 

/^yiWd. = 2.iog(7^) 


for /3 > 1. 
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Since it is expected that solutions to (2) converge to an equilibrium with the same 
mass if available, one anticipates that solutions with mass m < rric are globally 
dehned for all t G [0,cxo), while it is conjectured that solutions to (2) with mass 
m > nric will blow up at a certain hnite time T*. 

Despite its simple formulation, very little can be said about the Cauchy problem 
for the BEFP equation (2) beyond the local existence of solutions [8, 7, 33], which 
can be proved by usual hxed-point arguments. It was proved by Toscani [33], via 
a contradiction argument involving the second moment, that in dimension d > 3 
there are indeed certain solutions which blow up if the initial mass is sufficiently 
large or, for any given supercritical mass, if the initial second moment is small 
enough. The nature of the blow-up and its prohle is an open question in d > 3. 
The numerical approach becomes then of even greater interest in clarifying the long 
time asymptotics of quantum models related to Bose-Einstein distributions, see for 
instance [34, 26]. We hnally remark that in one dimension a direct minimization 
of the entropy functional for related equations can lead to singular parts in the 
equilibrium states [3]. 

We focus here on the 2-dimensional case, where the critical mass is inhnite 
and show that smooth solutions exist globally in time under certain mild conditions 
on the initial data. We always consider non-negative solutions to (2). One of the 
main difficulties that equation (2) presents when attempting to show the existence 
of global solutions is the lack of a priori estimates that would allow to control the 
L^-norm of solutions for any p > 1. 

Our analysis is based on the following change of variables, where, for the sake 
of clarity we shall omit, in general, the dependence on t and r or n of the different 
quantities that we will introduce. Consider a radially symmetric smooth solution 
f(t,v) = |n|) to equation (2), with hnite mass m = ||/o||li > 0. Then 

(f = (p(t,r) satishes, for t,r > 0, 

dt<p = dip + dr(^(^r- + dr . (4) 

We notice that, since / is bounded at v = 0, p satishes that limr^o r~’^p(t,r) = 0 
for all k < d — 1. Since / is radially symmetric and smooth, V/(f,0) = 0 and p 
also satishes the boundary condition 

dr {r^~‘^p) =0 at r = 0. 

Expanding the derivative one sees that in particular hm^_j,o drp{t, r) = 0. We now 
dehne $(f,r) as the indehnite integral of p: 

:= [ p{t,s)ds. (5) 

Jo 
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That is, dud^ is the mass of / in a ball of radius r, where Ud is the volume of the 
d-dimensional unit ball. One sees that $ satishes 

d — 1' 


dt^ = + 


r — 




( 6 ) 


<h(f, oo) 


m 

dud 


From now on, £x the dimension d = 2. Then the coefficient of the last term on 
the right hand side of equation (6) becomes independent of r and we can use the 
Hopf-Cole transformation. We introduce the new quantity 


= - 1 . 


( 7 ) 


Then, equation (6) becomes: 

{ (9t\h = + (r — 

V "7 (8) 

T(t,oo) = - 1, 

in which the quadratic term has been absorbed by the second derivative of T with 
respect to r. We finally introduce the change 


:= dr^{t,r), (9) 

which leads us to 

dt'ijj = + dr . ( 10 ) 

At this point we only need to notice that equation (10) coincides with the radial 
form of the 2D linear Fokker-Planck (FP) equation, namely 

dtg = Ag + div(n 5 f), f > 0, n G (H^-) 

g{0,v) = go{v), n G (11b) 


More specifically, if g is radially symmetric and satisfies (11), then the function 
^|J = 'ip{t,r) defined by = \v\g(t,v) solves (10). Hence, all the existing 

theory for equation (11) [29, 18] can be used and translated back to the BEFP 
equation (2) through the different changes of variable. This gives in particular an 
analytic expression for radially symmetric solutions of (2) and directly yields an 
existence and uniqueness theory for the (BEFP) equation (2) in dimension 2, as 
well as information about the asymptotic behavior of its solutions. 

The above changes of variable are summarized in the following diagram: 


/ 


rf=‘Pir) 




rr 

Jo 


<F 


'’-1 




dr'^ 




rg=^{r) 


f=m 
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The rest of the paper is organized as follows. In Section 2 we make precise the 
translation between solutions to the Bose-Einstein-Fokker-Planck equation (2) and 
solutions to the linear Fokker-Planck equation (11). Next, in Section 3 we present 
our main results on radially symmetric solutions, Theorems 3.3 and 3.6. Then, in 
Section 4 we extend these results to more general solutions by comparison with 
radial solutions, and state our main results for non-radially symmetric solutions 
in Theorems 4.2 and 4.3. We state the Csiszar-Kullback inequality for the Bose- 
Finstein-Fokker-Planck entropy in Lemma 4.5, which is needed for the proof of 
Theorem 4.3. 


2. The change of variables 


We will now make precise the equivalence between radially symmetric solutions 
of the FP equation (11) and those of the BFFP equation (2). This will hold 
for classical solutions of BFFP, as given in the next dehnition. We denote by 
the space of functions / = fit,v) on (0, cxo) x for which /, dtf, and all 
partial derivatives of order up to 2 with respect to v exist and are continuous on 
(0, cxo) X R^. A measure is always understood to be a non-negative Borel measure 
with hnite total mass and the space of such measures is denoted by A1+(R^). 

Definition 2.1 (Classical solution to BFFP). Let /o be a nonnegative Borel mea¬ 
sure in R^ (not necessarily hnite). We say that a function / dehned on (0, cxo) x R^ 
is a classical solution to the Cauchy problem for the BEFD equation (2) with initial 
condition /o when / = f{t,v) belongs to ClC^, equation (2a) holds in (0, cxo) x R^, 
and limt_,.o f{t) = /o in the weak-* sense of measures on R^. 

We dehne a classical solution to the Cauchy problem for the FP equation (11) 
analogously. Note that classical solutions of the FP equation are known to be 
unique for measures with hnite total mass as initial data via the change of variables 
to the heat equation, see next section and [9]. However, since our change of 
variables works for locally hnite measures, we prefer to state the concept in this 
more general setting. 

The change of variables (5), (7), and (9) can be summarized as follows: the 
non-negative, radially symmetric functions / and g on R^ are related by 


v\f{t,v) = ip{t,\v\), \v\g{t,v) t > 0, n e R^ (12) 


with 



for t > 0, r > 0, 


(13) 


or equivalently 


'ipit, r) := dr exp 



for t > 0, r > 0. 


(14) 
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This is also equivalent to 


ip{t,r) = 


1 + ds’ 


= ^p{t,r)exp i (^(t,s)ds , (15) 


which may be a more convenient form for some calculations. 


We may also adapt the change of variables (12)-(14) to work on locally hnite 
non-negative measures (i.e., measures which are hnite on any compact set), so that 
we can consider a measure as initial condition. If / and g are non-negative, locally 
hnite, radially symmetric measures on then one may dehne measures (p and ^|J 
on [0, cxo) by 

■= , ij{[0,r]) := , r > 0. (16) 

Conversely, if p and "0 are radially symmetric measures on [0, cxo), then there are 
unique radially symmetric measures / and g on R^ which satisfy (16). When / 
and g are functions, p and ^|J are just dehned by (12). 

The change (13)-(14) may be extended to non-negative, locally hnite measures 
by setting 

V9([0, r]) = log (1 -I- -^([0, r])) for r > 0, (17) 

-?/;([0, r]) = exp (g3([0, r])) — 1 for r > 0. (18) 


One can check that the change (16)-(18) coincides with (12)-(14) when / and g 
(and hence also p and ip) are locally integrable functions. 

Hence, by setting K{g) = / where / is related to g through p and j/^, see 
equations (16)-(18), we dehne a transformation A on the set of non-negative, 
radially symmetric, locally hnite measures on R^. We observe that from (17)- 
(18) one easily deduces that A and A“^ are continuous in the weak-* topology 
of measures. We will often consider A restricted to the set A4+^i.ad(lR^) of non¬ 
negative, radially symmetric, hnite measures on R^. In that set, the masses (total 
variation) of / (denoted by m) and g (denoted by M) are related by 

m = 271 log (^1 + , M = 271 _ i) . (ig) 

The transformation A is also continuous in Lj^j(R^), the set of integrable and 
radially symmetric functions in R^: 

Lemma 2.2. The change A is a bijection in L);jj^(R^), and both A and A~^ are 
locally Lipschitz continuous in the L^-norm. 

Proof. That A is a bijection in L):g^^(R^) is clearly seen from (15), for example. Take 
gi and g 2 in and set /* = A^gf), |n|v9i(|n|) = fi{v), and |n|V’i(|n|) = gi{v) 

for 7 = 1,2 and v G R^. Set Ti(r) = (1 + /J" 'ipiis) ds)“^ for i = 1, 2 and notice that 

0 < Tj(r) < 1 and \Fi{r) — F 2 {r)\ < f |'0i(s) —'02(s)| ds , r > 0. 

Jo 
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Combining these inequalities with (15) gives 

noo noo 

^ll/i - /2II1 = y - (P2{r)\dr = J \Fl{r)^jJl{r) - F2{r)i)2{r)\dr 

poo poo 

< / Fi{r) \'ijji{r) -'i/j 2 {r)\ dr + / \Fi{r) - F 2 {r)\'ijj 2 {r) dr 

Jo Jo 

poo poo pr 

< / |i/^i(r) - ■?/; 2 (r)| dr + / / |■^/;l(s) - ■?/; 2 (s)| i/^2(7’) ds dr 

Jo Jo Jo 

- ^ “ ^2||l||^2||l, 

which shows that A is locally Lipschitz continuous in the L^-norm. A very similar 
argument shows that A“^ is also locally Lipschitz continuous. □ 

This change of variables is also very regular, in the sense that it preserves the 
smoothness of the function to which it is applied. This will enable us later to 
translate properties of the FP equation (11) to the BEFP equation (2) in dimen¬ 
sion 2: 

Lemma 2.3. The change A is a bijection from the set of radially symmetric func¬ 
tions in C^(M^) to itself (in the sense that A{g) coincides a.e. with a function in 
C^(M ^)). In a similar way, if g E C(C^ with g(t) radially symmetric for each t > 0 
and we define f{t,v) = A{g(t)){v) then f G 

Before proving this we mention an elementary lemma that we wish to make 
explicit. We denote by C^([0, 00 )) the space of twice differentiable real functions 
on [0, 00 ) such that all derivatives of order up to 2 have limits at 0. 

Lemma 2.4 (Regularity of a radial function). Take 4> G F^([0, 00 )) such that 

lim (f'ir) = 0. (20) 

If we define a radial function f: ^ R by f{v) = 4>{\v\), then f G C^(M'^). 

Conversely, if f E C^(M'^) is a radially symmetric function, then the function 
<p\ [0, 00 ) —M defined by 0(r) := /(r, 0,..., 0) belongs to C^([0, 00 )) and satisfies 

( 20 ). 

Remark 2.5. It is easy to extend this lemma also to time-dependent functions: for 
example, if a function cf = (f>(t,r) is in C/C^((0, 00 ) x [0, cx))) and satishes 

lim dr(t>(s,r)=0 
(s,r)— 

for all f > 0, then f{t,v) := (f{t, |u|) is in ClC^. 

Proof of Lemma 2.f. Clearly / G C^(M'^ \ {0}) and we only need to check that the 
hrst and second-order derivatives of / can be extended continuously at u = 0. For 



8 JOSE A. CANIZO, JOSE A. CARRILLO, PHILIPPE LAURENgOT, AND JESUS ROSADO 


1 < *, j < and V G v ^ 0, 

dif{v) = > didjf{v) = 5ij(t)"{\v\) + (^0"(|n|) - • 

Since v i—)■ Vi/\v\ and v i—)■ ViVj/\v\'^ are bounded and 

lim (ct)"{r) - = 0"(O) - 0"(O) = 0 

r->0 V ^ / 

by (20) we conclude that 

\im. dif{v) = 0 and \im didjf{v) = (j)''{0)6ij . 

tl —>-0 D —>-0 

Consequently, f can be extended to a C^-smooth function in by setting V f(0) = 
0 and £>V(0) = 0"(O)id. 

The converse is obvious as the C^-regularity of 0 readily follows from that of / 
while the radial symmetry of / guarantees that V/(0) = 0 and thus (20). □ 


Proof of Lemma 2.3. Take a radially symmetric function g G C^(R^) and let us 
show that f = A{g) E C^(R^) (the proof for is similar). Clearly / is 
smooth in R^ \ {0}, since the change of variables has a smooth expression away 
from n = 0. In order to see that / is twice differentiable at n = 0, we use (12) and 
(15) to write / as 


fi'v) 


Ini 


V’(k]) 

|n| 



"hi 


^l){s)ds 


-1 




where we have set 


F(r) 



tjj{s)ds 


-1 


r > 0. 


Observe that '0(r) = rg{r, 0,..., 0) for r > 0 and thus belongs to C^([0, cxo)) with 
■0(0) = 0. then T'(0) = 0 and we deduce from Lemma 2.4 that / extends to a 
C^-smooth function in R^. The time-dependent results can be proved with very 
little modihcations. □ 


We have the following result; 

Theorem 2.6 (Equivalence between BEEP and FP). Let fo and go be non-negative 
radially symmetric measures on R^ related by (16)-(18), and let f, g be non¬ 
negative radially symmetric functions on (0, cxo) x R^ related by (12)-(14). The 
following are equivalent: 

(1) f is a classical solution to the Bose-Einstein-Fokker-Planck equation (2) 
with initial condition fo. 

(2) g is a classical solution to the linear Fokker-Planck equation (11) with 
initial condition go. 
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Proof. We just show one of the implications, since the reverse one is analogous. 
Assume point ( 2 ), and let us show point (1). Since 5 ^ is a classical solution to (11), 
it lies in Hence the calculations (4)-(10) are rigorous and we see that / 

satishes ( 2 ) in ( 0 , 00 ) x (M^ \ { 0 }) (the point u = 0 being a possible problem due 
to the radial change of variables.) But due to Lemma 2.3, / belongs to and a 
continuity argument entails that / satishes (2a) in (0, 00 ) x Finally, / satishes 
the initial condition in the weak-* sense since the change A is continuous for the 
weak-* convergence of measures. □ 

Notice that some of the classical solutions given by Theorem 2.6 are not covered 
by any of the existence results in the literature. The change of variables gives 
us classical solutions to a nonlinear equation with a measure as initial condition, 
which is quite remarkable. Finally, let us point out that radially symmetric classi¬ 
cal solutions to the FP equation (11) are unique as soon as some growth condition 
at inhnity is satished. Since the FP equation (11) is equivalent to the heat equation 
as we will see next, then uniqueness is satished for initial data in tempered dis¬ 
tributions (classical solutions are distributional solutions), in particular for hnite 
mass non-negative measures (which are our main interest). 

3. Global radially symmetric solutions 

The Fokker-Planck equation (11) is equivalent to the heat equation via a self¬ 
similar change of variables. Moreover, it can be explicitly solved through the 
Fourier transform and its fundamental solution can be explicitly found leading to 

J^{t,v,w) := a(f)"'^/^M^(q(a(f)"^/^u - w) 

with 

a{t) := e- 2 * , io{t) := e'* - 1 and Ma (0 := ( 27 rA)-''/ 2 g-|qV 2 A 

for any A > 0. Let us focus again on the two dimensional case and dehne the 
operator acting on functions f by: 

J^[.^(tn)](f, u) := / T{t,v,w)f{w) dw , (f, u) G ( 0 , cxo) x . ( 21 ) 

The following result is well-known, see for example [9, Lemma 3.1]: 

Theorem 3.1 (Solutions to FP). Let ^'(R^) denote the space of tempered dis¬ 
tributions. Then, given go G S"(R^) there exists a unigue weak solution g G 
C([0, cxo), S"(R^)) to the linear Fokker-Planck eguation (11) with initial condition 
go given by g = J^[go\- In addition, g is smooth on ( 0 , 00 ) x R^, and positivity and 
radial symmetry are conserved; more precisely: 

(1) If go is not identically egual to zero then g{t,v) > 0 for all t > 0 and 
V G R2. 

(2) If go is radially symmetric then g{t) is radially symmetric for all t > 0 . 
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Remark 3.2. Since endowed with the weak-* topology is continuously 

embedded in in our context this result tells us that, indeed, for any qq e 

there exists a unique classical solution g G C([0, cxo), of the FP 

equation (11) (where we always consider the weak-* topology on A1+(M^)) and 
this solution satishes the above properties. 


Finally, let us recall some moment and smoothing properties of the solution of 
the FP equation (11). With this aim, let us introduce the spaces ^^(M^) with 
norms 

\\f\\q-=\\{^ + \v\^)f\\p and \\f\\p:= (^J^Jf\Pdv^ . 

for 1 < p < oo and i > 0, with the usual essential supremum norm dehnition for 
the case p = oo. Solutions g to equation (11) satisfy 


\\dag{t)\\q < 




for all 1 < g < p < oo, £ > 0, a G a multiindex, and t > 0. Here, C* > 0 is 
a number depending on p, q, £ and a. For a proof of this result see [7, Proposi¬ 
tion A. 1]. In particular, the case a = 0 and p = q gives 

\\gmq<Ce'-^\\go\\q (22) 

for all 1 < p < oo, £ > 0, and t > 0 with C depending only on p and £. The case 
a = 0, g = 1 gives 

\m\\L;<Cq^] IlgolL; . «>0, (23) 


Theorem 2.6 immediately implies an analogous result for radially symmetric 
solutions to the BEFP equation (2): 

Theorem 3.3 (Radial solutions to BEFP). For any radially symmetric /o G 
A1+(]R^) there is a unique classical solution f G C([0, cx)), A1 + (M^)) to the BEFP 
equation (2) with initial condition fo and it is given by f{t) = A{g{t)) for allt > 0 
where g is the solution to the FP equation (11) with initial condition go := dV~^{fo) ■ 
In addition, f is bounded in L°°{{5,oo) x M^) for each <5 > 0 and fit) is radially 
symmetric for all t > 0. Furthermore, if fo is not identically equal to zero then 
f{t, n) > 0 for allt > d and n G M^. 

Remark 3.4. This solution is also the unique mild solution in the sense of [7, 
Eq. (2.2)] for initial data in the class considered therein, since the authors prove 
that mild solutions are unique up to the maximal existence time. Of course, the 
above theorem also shows that the unique mild solution does not blow up and is 
dehned for all positive times. 
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Proof. Consider the solution g to the FP equation (11) with initial condition 
A“^(/o) given by Theorem 3.1. This solution is also in C([0, cxo), (see 

Remark 3.2). Since is a classical solution to the FP equation (11), the function 
/ ;= A{g) is a classical solution to the BEFP equation (2) by Theorem 2.6. It 
also lies in C([0, oo), A1+(M^)), since A preserves this space and is continuous with 
respect to the weak-* topology of measures. Finally, the claimed boundedness in 
(5, cxd) X readily follows from (23). □ 

We next give another example of the results that one can get by translating 
known results for the FP equation (11), which deals with the large time behavior. 
The following result is well known and can be found for example in [9, Theo¬ 
rem 4.1]. 

Theorem 3.5 (Asymptotic behavior for FP). Take go G L],_(M^) and denote the 
total mass of go by M. The solution g of the FP equation (11) with initial condition 
go satisfies 

|| 5 f(f) - MfifoolUi < e~^\\go - Mgc^Wii for all t > 0, 
where the equilibrium g^o is given by 

9ooiv) = V G 

ZTT 

Notice that if we apply the change (12)-(13) to Mg^o we obtain precisely /^, 
with 

,3 = 1 + 1 . ( 24 ) 

Thanks to Theorem 2.6 and Lemma 2.2, this translates into the following result 
for the BEFP equation (2) 

Theorem 3.6 (Asymptotic behavior for radial BEFP). Take a radially symmetric 
fo G L]^(M^) and denote its mass by m. There is a constant K = K{m) > 0, 
depending only on m, such that the solution f of the BEFP equation (2) with 
initial condition fo satisfies 

ll/(«) - < A+-‘ll/o - ftL\\L> for oil t > 0, 

choosing (3 to be such that f^ has mass m: 

/ 3 = (25) 

Remark 3.7. The constant K{m) above can be taken to be 

K{m) = (Lipjv^(A)) (Lip^(A-^)) , 

where Lip^ denotes the L^-Lipschitz constant in the L^-ball of radius R (see 
Lemma 2.2.) 
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Another interesting consequence of Theorem 2.6 is that radially symmetric so¬ 
lutions to the BEFP equation (2) have an explicit expression in terms of the initial 
condition /q, obtained by transforming the explicitly computable solution g of the 
FP equation (11) through the change of variables (12)-(15): 






(f, v) e (0, cxo) X 


1 -h /j’"' s) ds 

with ^|J related to g by (12). This relation obviously implies that, for {t,v) G 
(0, cx)) X 

271 


2vr -h IIs'oIIli 


g{t,v) < f{t,v) < g{t,v) < 


by the conservation of the L^-norm. This readily translates bounds for g into 
bounds for /; indeed, on the one hand, equation (22) gives 


||/(t)||Lf < Ce ’ II^oIIlp < Ce ^ UoUlp (26) 

for all 1 < p < oo, £ > 0, and f > 0 with C depending only on p and £. On the 
other hand, equation (23) gives 


ll/(t)llL? 


< 



< c 



ip-l)/p 

IIs'oIIli 


(p-i)/p 

gll/oIhlATT 



(27) 


These bounds on / will be used in the next section in order to show that solutions 
to (2) starting from general inital data are globally dehned in time. 


An interesting particular solution is the transform of the solution to the FP 
equation (11) with initial condition the Dirac mass ho at n = 0 which may be 
referred to as a fundamental solution for the (nonlinear) BFFP equation (2): 

1 r 1-1 


/(<■ *’) = 5(0 -1 


for V G M , t > 0, 


with 

dit) ;= (1 — e“^*) for t > 0. 

There is also an explicit inhnite-mass solution to the BEFP equation (2) ob¬ 
tained by transforming the solution g{t,v) = A > 0, to the FP equation 

( 11 ): 

f(t, v) = 2 (2A“^e“^* -I- |n|^) ^, for n G t > 0. 

Strictly speaking, this solution does not satisfy Dehnition 2.1 because it is not 
integrable; however, it satishes (2a) for all {t,v) G (0, cx) x This solution 
converges exponentially fast to an inhnite-mass stationary solution f^{v) := \v\~‘^. 
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4. Global (non-radially symmetric) solutions 

In the case of Fermi-Dirac interactions, where the nonlinear drift in (2) is neg¬ 
ative, a local well-posedness theory and L^-contraction and maximum principles 
were proved in [7]. However, the proofs are totally independent of the sign of the 
drift and can be directly translated to our case: 

Proposition 4.1 (Summary of results from [7]). Let fo G {L^ fl be non¬ 

negative, with p > 2 and £ > 1. Then: 

Local Existence: There exists a maximal time T G (0, oo] depending only 
II/oIIl^ II/o||li such that there is a unigue solution f mC([0, T), (L^n 
L^)(M^)) of the BEFP eguation (2) which is non-negative and satisfies 

ll/(^)IUi = II/oIIli ) t G [0,T) . 

IfT < oo then ||/(t)||LP —> oo as t ^ T. 

Contractivity and Comparison Principle: Let T G (0, oo] and fi, /2 G 

C([0,T), (L^ n L^)(M^)) be two solutions to the BEFP eguation (2) with 
initial data fi^ and / 2 ,o in nL^)(R^), respectively. Then 

— /2(^)||li < Wfo ~ 5'oIIli , ^ ^ [0, T) 

Furthermore, if fip < / 2 ,o then fi(t,v) < f 2 (t,v) for all t G [0,T) and 
u G 


Thus one can deduce the following result from Theorem 3.3 and Proposition 4.1: 


Theorem 4.2. Take a non-negative initial condition fo G with i > 2. 

Then there is a unigue solution f G C([0, oo), L^(M^)) (in the sense of [7]j of the 
BEFP eguation (2) with initial condition fo. In addition, if\\fo\\L\ < oo, there is 
some C > 0 which depends only on fo such that 


II/WIIli + ll/WlUr < C for all t > 1. (28) 


Proof. The local well-posedness is provided by Proposition 4.1, so we only need 
to show that the norm ||/(t)|lL|° does not blow up in hnite time. Let T be the 
maximal existence time of / given by Proposition 4.1. We consider the radially 
symmetric solution ftt to (2) corresponding to the initial condition 



ll/olUf 

1 -i- \vY ’ 


V eMf , 


which we know to be globally dehned by Theorem 3.3, since f^ has hnite mass 
due to the condition i > 2. Clearly, fo < fo and by the comparison principle 
in Proposition 4.1 we know that f{t,v) < fttft^v) for all t G [0,T) and v G M^. 
Equation (26) with p = oo gives us 




#1 



14 JOSE A. CANIZO, JOSE A. CARRILLO, PHILIPPE LAURENgOT, AND JESUS ROSADO 

for some C > 0 which depends only on I and ||/o||li- This estimate prevents the 
L^-norm of /(t) to explode in finite time, so the solntion mnst in fact be globally 
dehned, that is T = oo. In order to obtain the additional statement (28) it snffices 
to apply eqnation (27) with p = 1 and p = oo. □ 


We now tnrn to the large time behavior of solntions to (2) and recall that 
eqnation (2), in any dimension d, has an entropy fnnctional given by 


H{f) = 


-f{v) + f{v) log/(n) - (/(n) + 1) log(/(n) + 1) 


dn 


(29) 


in the following sense: any solntion / to (2) with a bonnded moment of order i > 2 
satisfies 

2 


df 


H{f ) = -D{f) --//(! + /)(«’) t> + V log (j^) W 


dn. 


(30) 


The availability of this fnnctional leads ns to the following resnlt: 


Theorem 4.3. Take a non-negative /o G L“(M^) fl L^(M^) with i > 2. Then the 
solution f of the BEFP equation (2) with initial condition /o given by Theorem 4-2 
satisfies 

lim \\f{t) - = 0 , 


1—>-oo 


where fi is given by (25) in terms of the mass m of fo. 

Remark 4.4. Dne to Theorem 4.2 solntions to (2) emanating from non-negative 
initial data in L“(M^) fl Lj(M^) with i > 2 are bonnded in L°° globally in time 
and a simple interpolation argnment shows that for all 1 < p < cxo we also have 

lim \\f{t)-f^\\p = 0 . 

r—>-oo 


Proof of Theorem 4-3 ■ We follow a common techniqne to dednce convergence to 
eqnilibrinm from the existence of a Lyapnnov fnnctional snch as (29); see for 
example [11, 4]. 

Take a solntion / = f{t,v) of the BEFP eqnation (2) with initial condition 
/o and recall that Theorem 4.2 gnarantees that {f{t) : t > 1} is bonnded in 

Lf (M^) nLi(M2). 

Step 1: Strong convergence of f(tn) for some sequence {tn}n>i of times. 
From (30) (which we can nse since we are assnming /q G //^(M^) for i > 2) and 
the fact that H{f{t)) is bonnded below by H{f^) (dne to Lemma 4.5 below) 
we see that there mnst exist a divergent seqnence {tn}n>i, > 1, snch that 
-^(/(^n)) ^ 0 as n —)■ cx). Dne to (28) we know that there exist a snbseqnence 
of {tn}n>i (which we still call {t„} by an abnse of notation) and a non-negative 
fnnction f^ G fl snch that f{tn) foo weakly in L^(M^), for all 
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p G [l,cxo), and also in the weak-* topology of We also note that the 

bounds (28) ensure that 


/oo(n)dn = lim / /(t„, u) dn = m := ||/o||li 


(31) 


Setting fn := f(tn), n > 1, let us now see that the convergence of {fn}n>i must 
in fact be strong, due to the fact that D{fn) —?■ 0. Indeed, expanding the square 

r r |v/„i2 


Difn) = fnil + fn)iv)\v\‘^ dv + 

2r2 

IV/J2 


/n(l + fn) 


{v)dv + 2 V ■ Vfn{v) dn 


> 


/n(l + fn) 


(v) dn — 4m. 


Since D{fn) —)■ 0, this shows that the sequence {log(/„/ (l + /n))}n>i is bounded in 
H^(Mf). In particular, there is a subsequence of {fn}n>i (again denoted by {fn}) 
which converges pointwise a.e. to /oo and weak limits coincide. Vitali’s theorem 
then shows that, after possibly extracting a further subsequence, we can conclude 
that 

fn foe strongly in L^(Ef) as n —)> oo, 

for all 1 < p < oo. 

Step 2: Proof that f^o = f^- Dehne hn '■= /n/(l + fn), n > 1. From the 
previous step, we know that hn converges pointwise a.e. to hoc ■= /oo/(l + /oo)- 
Since hn < fn is clearly uniformly bounded in all L^iMf) {p G [l,oo]), there is 
another subsequence of {fn}n>i such that {hn}n>i converges weakly in all L^iMf), 
and weakly-* in L°°{Mf) to h^. Now we have, by the Cauchy-Schwarz inequality. 


\vhn{v)+ Vhn{v)\ dv] < 



V -F 




II II Ll 


< mD{fn) -t 0 


as n —)■ oo. This implies that vhn + Vh„ —?■ 0 in L^(M^) as n —?■ oo which, to¬ 
gether with the already established convergences of hn to hoo implies that vhoo{v) + 
Vhoo(u) = 0 in the sense of distributions and in fact almost everywhere. Conse¬ 
quently, hooiv) = /3exp(—|np/2) for a.e. n G for some /3 > 0 and thus /oo = /^ 
a.e. in R^. Finally, f3 = f3 is given by (25) since the mass of /oo is the same as that 
of /o by (31). 

Step 3: Convergence of the entropy. Since t e-)■ H{f{t)) is nonincreasing and 
larger than H{f^) (see Lemma 4.5 below) there is Hoo > H{f^) such that 

H{f{t)) —> Hoo as t —)■ oo. (32) 

From the previous two steps we know that {/„ = /(tn)}n>i converges pointwisely 
to /^ while (28) ensures that {fn}n>i is bounded in {Lj fl L^)(R^), hence also 
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in Since £> 2 the Dominated Convergence Theorem now shows that 

H{fn) -)■ which implies that 

Using Lemma 4.5 we directly dednce from (32) that ||/(t) — /^||i ^ 0 as t —?■ oo, 
hnishing the proof. □ 


We hnally give a Csiszar-Knllback ineqnality for our entropy functional, which 
was needed to show that convergence in entropy implies convergence in L} norm in 
the proof of Theorem 4.3. We would like to warmly thank the anonymous referee 
for suggesting an improvement of the following result in the preprint version of 
this paper. 

Lemma 4.5 (Csiszar-Kullback-type inequality for the BEFP entropy). Take a 
non-negative function f G L^(M'^) with mass ||/||i < me, and (3 > 0 such that the 
mass of is equal to that of f. Then 

(33) 

for some C > 0 which depends only on ||/,^||i and ||/^|| 2 - 


Proof. We follow a well-known strategy for proving Csiszar-Kullback-type inequal¬ 
ities; see for example [5]. For each n G call 

InP 

(j)(r) = (j){r, v) := r log r — (r -|- 1) log(r -|- 1) -|- r > 0, 

so that we have 

H{f)= [ Hf{v))dv. 
jRd- 

Using that 4>'{f^{v)) is a constant independent of n G we can write 


H{f) - 

= (0(/(w)) - HfLio)) - (p'{fL{v)){f{v) - /^(u))) dv. (34) 

Now we notice that 0 is a convex function with 4>''{r) = (r(l -|- n))“^, so we have 
the bound 

0(r) - 0(s) - 0'(s)(r -s)> > 0 

s(l + s) 

whenever 0 < r < s. Using this in (34) and denoting D := {n G | f{v) < f^{v)} 
we obtain (writing for convenience / = f{v) and f^ = f^{v) inside the integral) 




> 


if - fLr 

foo{l + foo) 


dn > / (/^ - /) dn 


fL{o){l +f^)dv 


-2 
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where for the last step we used the Cauchy-Schwarz inequality. Calling 

and observing that 

f (fL-f)'iv = h\f^-f\U 

tj 

this gives precisely 

□ 
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